Monte Carlo study of an unconventional superconducting phase in Ir-oxide Jeff = 1/2 

Mott insulators induced by carrier doping 
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Based on a microscopic theoretical study, we show that novel superconductivity is induced by 
carrier doping in layered perovskite Ir oxides where a strong spin-orbit coupling causes an effective 
total angular momentum Jeff = 1/2 Mott insulator. Using a variational Monte Carlo method, we 
find an unconventional superconducting state in the ground state phase diagram of a t2g three-orbital 
Hubbard model on the square lattice. This superconducting state is characterized by a d^2_y2-wsive 
"pseudospin singlet" formed by the JcS = 1/2 Kramers doublet, which thus contains inter-orbital 
as well as both singlet and triplet components of t2g electrons. The superconducting state is found 
stable only by electron doping, but not by hole doping, for the case of carrier doped Sr2lr04. We 
also study an effective single-orbital Hubbard model to discuss the similarities to high-Tc cuprate 
superconductors and the multi-orbital effects. 

PACS numbers: 74.20.-z, 71.30.-|-h, 75.25.Dk 



Search for novel superconductivity (SC) is one of the 
most interesting and fundamental issues in condensed 
matter physics. In strongly correlated electron systems, 
SC is very often induced in the vicinity of long-range 
ordered states, thus suggesting the importance of the en- 
hanced fluctuations for the SC. The most studied exam- 
ple is found in high-Tc cuprate superconductors [l| where 
the SC occurs next to an antiferromagnetic (AF) Mott 
insulator, and thus the AF spin fluctuations are often 
believed to be responsible for the SC 0- 

RecentlVjlayered perovskite 5c? transition metal oxides 
Sr2lr04 0,11 andBa2lr04 have attracted much atten- 
tion because several experiments have revealed a novel 
spin-orbit-induccd Jeff = 1/2 Mott insulating behavior 
at low temperatures @-[l3|- l^ these systems, due to a 
large spin-orbit coupling (SOC) and a large crystal field 
splitting, the local electronic state with nominally {t2g)^ 
electron configuration in Ir ion is represented by an effec- 
tive total angular momentum Jcs^\—L + S\ = l/2 p^ . 
In these insulators with effectively one hole per Ir ion, this 
pseudospin remains a good quantum number and orders 
antiferromagnetically. Indeed, very recent experiments in 
Sr2lr04 have observed that the low-energy magnetic ex- 
citations can be well described by a pseudospin 1/2 AF 
Heiscnber g rn odel with AF exchange coupling Jcx=60- 
100 meV (lol. [T^. The theoretical studies also support 
this Jeff = 1/2 Mott insulator in these Ir oxides [14l - [l8| . 

It is now interesting to compare Sr2lr04 with parent 
compunds of high-Tc cuprate superconductors such as 
La2Cu04. The similarities are summarized as follows: 
(i) both are in the same layered perovskite structure 
of K2NiF4 type, i.e., in a quasi two-dimensional (2D) 
structure, (ii) both have effectively one hole per Ir or Cu 
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ion, (iii) both show spin or pseudospin 1/2 AF order at 
low temperatures with low-energy magnetic excitations 
described by a spin or pseudospin 1/2 AF Heisenberg 
model, and (iv) both have the large Jex of the same or- 
der. Considering these similarities, it is tempting to ex- 
pect that novel unconventional SC with possibly a high 
critical temperature (Tc) is induced once mobile carri- 
ers are introduced into the Jeff = 1/2 Mott insulating 
Sr2lr04. Although there have been several reports sug- 
gesting this possibility [l3, El] , it is highly desirable to 
show, based on microscopic calculations, the existence of 
SC in the doped Jeff = 1/2 Mott insulator. 

In this Letter, using a variational Monte Carlo (VMC) 
method, we study the ground state phase diagram of a 
three-orbital Hubbard model, and show that novel 
unconventional SC is induced by carrier doping in the 
JcS = 1/2 Mott insulator. This SC is characterized 
by a dj;2_y2-waYe "pseudospin singlet" formed by the 
Jc« = 1/2 Kramers doublet, which thus consists of inter- 
orbital pairings and both singlet and triplet pairings of 
t2g electrons. We also find that the SC is stable only by 
electron doping, but not by hole doping, for the case of 
carrier doped Sr2lr04. We furthermore study an effective 
single-orbital Hubbard model to discuss the similarities 
to high-Tc cuprates and the multi-orbital effects. 

One of the simplest models for the Ir oxides, which we 
shall study here, is a t2g three-orbital Hubbard model 
on the square lattice defined hy H = Jfkin + ^so + Hi, 
where Jfkin = J2k,a,a ^ai^)^iaa(^kaa IS the kinctic 
term, i?so = X^i -^i ' the SOC term with a 

coupling constant A, and Hi = ?7 a "■icT'^iaJ- ^" 

^i,a<p,a [U'riiaaTlipg + ([/' - J)?!^^^"*/?^] + 
JY.i,a<p{c\^Afii'^^°'^'^^P^ + 4aTcL;Ci/3iCi^it + H.C.) 

is the Coulomb interaction term including intra-orbital 
([/), inter-orbital [U'), and spin-flip and pair-hopping 
(J) interactions [l^. Here, c\^^ is an electron ere- 
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FIG. 1: (a) Fermi surface and (b) energy dispersions of 
the non-interacting tight-binding energy band for Sr2lr04 
with electron density n — 5. Numbers in (b) denote 
the band index m, and Ef is the Fermi energy. A set 
of tight-binding parameters used is {t2,t3,t4,t5, jj.xy, X) — 
(0.5, 0.25, 1.03, 0.17, -1.0, 1.39)ti. 



ation operator at site i with spin a{=^,],) and orbital 
a(= yz,zx,xy), riiaa = c\a„c.iaa, and c\^^ is the Fourier 
transformation of c] . We impose U = [/' + 2J for 
rotational symmetry [20|. 

The kinetic and the SOC terms can be combined, 
Ho(ti, jixy, A) = -ffkin + ^^sOj in the matrix form 



kyza ' kzxa ' kxya 



£yz{k) iScrA/2 -S„\/2\ / Ckyza 

X I -is^A/2 e^:j(fe) iA/2 Ckzxo 
-.ScrA/2 -iA/2 exy{k) J \ckxyg. 



(1) 



E 

k.m.s 



where a is the opposite spin of a and = 1(— 1) 
for a =t (i). Notice that the SOC mixes the differ- 
ent electron spins (a and a), and the new quasipar- 
ticles, obtained by diagonalizing Hq, are characterized 
by band index m{= 1,2,3) and pscudospin s = (ti-l) 
with a creation operator a^^^. In the atomic limit with 
£yz{k) = Szx{k) = £xy{k) = 0, the sixfold degenerate 
t2g levels are split into twofold degenerate J^s = 1/2 
states (m = 1) and fourfold degenerate J^s = 3/2 states 
(m = 2, 3) [13. The undoped filling corresponds to elec- 
tron density n = 5, and in the atomic limit all states but 
the Jeff = 1/2 states are fully occupied. 

In Ref. Ilg . we have constructed the non-interacting 
tight-binding energy band for Sr2lr04: £yz{k) = 
—2^5 cos kx — 2^4 cos ky, ezx{k) = —2^4 cos k^ — 2^5 cos ky, 
and Exyik) — ~2ti{coskx +cosky) — 4^2 cos fc^ cos fcj, — 
2^3 (cos 2/02; + cos2ky) + fixy with a set of tight-binding 
parameters (^1,^2, ^3, ^4, t^^fixy, A)=(0.36, 0.18, 0.09, 0.37 
0.06, —0.36, 0.50) eV. The corresponding Fermi surface 
and energy dispersions are shown in Fig. [TJ As shown 
in Fig. [1] (b), we assign the band index m = 1, 2, and 3 
from the highest band to the lowest one, and only band 
1 crosses the Fermi Energy [U . 

The effect of Coulomb interactions is treated using a 



VMC method [Ig]. The trial wave function l^") consid- 
ered here is composed of three parts: j^*) = Pj^Pq |^). 
1$) is the one-body part obtained by diagonalizing Hq = 
Holti, flxy, Xap) with variational "renormalized" tight- 
binding parameters {ti, flxy, Xap}- Notice that we in- 
troduce an orbital dependent "effective" SOC constant: 
A Aq^. To treat magnetically ordered states, a term 
with a different magnetic order parameter is added to 
Hq. Here, we consider out-of-plane AF order (along z 
axis, 2-AF) and in-plane AF order (along x axis, x-AF), 
described respectively by I]j,m ^me"^ ''' (al„-,.aimt - 

is the Fourier transformation of a^.^^ and 



where a, 

Q = (tTjTt). The order parameters (Aff, M|,M|) for 
z-AF and (Aff , Mf, Mf) for x-AF are variational pa- 
rameters. With an appropriate basis transformation, we 
obtain the original t2g orbital representation in real space 
and construct the Slater determinant |$) for VMC sim- 
ulation. 

To study a possible superconducting state, we consider 
the following BCS-type Hamiltonian, 



-^BCS — E ("fclT' ^fc2t' °^fc3t' "-fcU' 0.-k2l, a-k3lj ^ 
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where £,m = Em{k) — /i, Em{k) is the eigenvalues of Hq, 
and /i is a variational parameter of the chemical poten- 
tial form. The gap functions Amm' [k dependence im- 
plicitly assumed) are additional variational parameters. 
After diagonalizing Hbcs- the ground state of the su- 
perconducting state is obtained by creating all negative 
energy states (7J) and annihilating all positive energy 
states (7^) on the vacuum state, |<I>) = Hi/ 7i'7!'' 1^ 
this study, we consider mostly the intra-band (but inter- 
orbital) pairing, namely, A„i,n' = for m ^ m! p3| . 
(3) 

The operator Pq ' is a Gutzwiller factor extended for 
the three-orbital system and the operator Pj^ is a long- 
rang charge Jastrow factor. These operators are exactly 
the same ones reported in Ref. [ll]. The ground state 
energies are calculated with a VMC method. The varia- 
tional parameters, as many as 80 parameters for a 20x20 
square lattice, are simultaneously optimized to minimize 
the variational energy by using the stochastic reconfigu- 
ration method [23| . 

Let us first summarize in Fig. [5] our main results 
for the ground state phase diagram where the electron 
density n and the intra-orbital Coulomb interaction V 
are varied. In this phase diagram, we set the Hund's 
coupling J = 0. Therefore, the pscudospin rotational 
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FIG. 2: (color online) The ground state phase diagram of the 
three-orbital Hubbard model on the 2D 20x20 square lattice 
with J = 0. PM, AFI, AFM, and SC denote paramagnetic 
metal, AF insulator, AF metal, and d^2 _y2-weLYe "pseudospin 
singlet" SC, respectively. 



symmetry is preserved [l^ [l6j , and z- AF and x-AF arc 
energetically degenerate (denoted simply by AF in the 
phase diagram). At n = 5, the AF insulator appears 
for U/ti > 2.6. This insulating state is considered to 
be the spin-orbit- induced Jeff = 1/2 Mott insulator 
observed in Sr2lr04 [l^. With increasing n by electron 
doping, the AF order is eventually destroyed, and it is 
replaced by the superconducting state for U/ti > 6 and 
n ~ 5.2 (namely, ^ 20% electron doping). The Cooper 
pair symmetry of this SC is found to be a dx2-y2-wave 
"pseudospin singlet" formed by the Jeff = 1/2 Kramers 
doublet |25|. Thus, this pairing contains inter-orbital 
components as well as both singlet and triplet com- 
ponents of t2g electrons. This can be easily seen in 
the limit of large SOC (A — > oo), where the most 
dominant Cooper pair is expressed by ifcij-aLfci^ oc 



(c^ 



kyzX 



(c^ 

y -kxyi 



-kyzt 



+ ic' 



kzxf 



When A is finite, the coefficient of each term is different 
from the one in the above limit and it is determined 
variationally. It should be also noted that the supercon- 
ducting order parameters considered in Eq. ([2]) include 
not only the nearest-neighbor pairing but also long-range 
pairings, e. g., u p to the 5th neighbor for the d^2_y2--wavc 
symmetry [26|. It is known that the simplest dx^-y^- 
wave pairing (oc cosfc^^ — cos ky) is significantly modified 
by the long-range contribution in the underdoped regime 
of high-Tc euprates (27j . The long-range pairings are 
important also in our three-orbital model and give the 
lower variational energy for the superconducting state. 

To the contrary, with decreasing electron density by 
hole doping for n < 5, we do not find a superconduct- 
ing state in the phase diagram (not shown here), where 
instead a AF metal dominates the SC [l^. This electron- 
hole asymmetry of the phase diagram reminds us of the 
phase diagram of a model for high-Tc euprates. As in the 
case of high-Tc euprates, the asymmetry found here for 
the three-orbital Hubbard model is understood due to a 
band structure effect. As seen in Fig. [T] (b) (indicated 



FIG. 3: (color online) The ground state phase diagram of the 
three-orbital Hubbard model on the 2D 20x20 square lat- 
tice. The solid-red and dotted-blue lines represent the phase 
boundaries for U/ti = 8 and 10, respectively. PM, a;- AFI, 
x-AFM, and SC denote paramagnetic metal, in-plane AF in- 
sulator, in-plane AF metal, and d^2 _y2-wa,ve "pseudospin sin- 
glet" SC, respectively. 



by the dotted circle), the dispersion of band to = 1 is 
fiat for k around (±7r, 0) and (0,±7r), which induces van 
Hove singularity in the density of states (DOS) for n > 5 
(but not for n < 5). The large DOS originated from this 
flat dispersion favors the SC [1^ . 

Next, we study the effect of Hund's coupling J. Fig. [3] 
shows the ground state phase diagram in a n-J/U plane 
for U/ti = 8. The introduction of a finite J breaks the 
pseudospin rotational sysmetry [13, [ill , and the states 
with in-plane AF order are favored over those with out- 
of-plane AF order (see Fig. [3]). We also find in Fig. [3] 
that the superconducting state remains stable, but the 
superconducting region gradually reduces and eventually 
disappears with further increasing J/U. This implies 
that the Hund's coupling J unfavors the SC. To under- 
stand the effect of Hund's coupling J, we recall that the 
charge gap Ac in the limit of strong Coulomb interactions 
with d^ configuration is Ac = E{d'^) + E{d^) - 2E{d^) = 
U — 3J [1^. Thus, the effect of the Hund's coupling J 
is to reduce the effective electron correlations. Compar- 
ing the variational energies of the paramagnetic metal 
and the superconducting state, we find that the SC is 
stabilized by the gain of the interaction energies at the 
expense of the band energies. Therefore, with increasing 
the Hund's coupling J, the condensation energy of the 
SC is greatly reduced and the SC is eventually destabi- 
lized [3l[. We also study the J/U dependence of the SC 
for U /ti = 10, and the results are shown in Fig. |31 Al- 
though the Hund's coupling J is still destructive for the 
SC, the superconducting region is found to be rather ex- 
tended as compared with the results for U /ti = 8. This 
suggests that the SC is more likely to be found experi- 
mentally in the doped Ir oxides with larger U/ti. 

It should be noted here that the realistic values of 
the Coulomb interactions for the Ir oxides are still con- 
troversial, and that even these values may vary greatly 
for different models because the screening mechanisms 
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FIG. 4: Energy comparison of the effective single-orbital Hub- 
bard model on the 2D 14x14 square lattice with U/t = 13, 
t'/t = 0.26, and t" /t = -0.05. The energies for different 
states are measured from the one for the paramagnetic state. 
The phase diagram is indicated in the upper part of the figure. 
Here n = 1 corresponds to half filling. PM, AFI, AFM, and 
SC stand for paramagnetic metal, AF insulator, AF metal, 
and d^2_y2-waYe SC, respectively. 

can be different [11] • Nevertherlcss, the previous stud- 
ies [13, m, Hi] have reported that U = 2-3 eV (i.e., 
U/ti = 5.6-8.3) and J/U = 0.05-0.20, which arc still 
within (or at least in the vicinity of) the parameter re- 
gion where we find the SC induced by electron doping. 

Finally, let us study an effective single-orbital Hub- 
bard model [m, [l^ to discuss the similarities to high- 
Tc cuprates and the multi-orbital effects. The effec- 
tive single-orbital model is readily constructed by fit- 
ting band m = 1, i.e., "Jeff = 1/2 band", in Fig. [T] 
(b) using the dispersion relation e{k) = —2t{coskx + 
cosky) — At' cos kx cos ky — 2t"{cos2kx + cos2ky) with 
{t,t',t") = (0.221,0.057,-0.011) eV. The value of U/t 
for the effective single-orbital Hubbard model is chosen 
to be 13 (w 8ti/t), which should correspond to U/ti =8 
for the three-orbital Hubbard model. 

The ground state phase diagram of this effective 
single-orbital Hubbard model is studied using the VMC 
method. The trial wave function used is almost the same 
as the one for the three-orbital model: j^*) = Pj^Pq"^ |$). 

The operator Pq"^ = Yii [1 ~ (1 ~ 5)'^it"-ii] ^ usual 
Gutzwiller factor with only one variational parameter g. 
For the one-body part |$), we consider the BCS-typc 
wave function with both AF and c?j.2_,y2-wavc supercon- 
ducting orders [35j . 

Fig. m shows the doping dependence of the variational 
energies for different states in the single-orbital Hubbard 
model. At half filling (71 = 1), the ground state is an 
AF insulator. By electron doping with n > 1, the AF or- 



der vanishes around 12% doping and the dx2-y2-wave SC 
appears. To the contrary, the SC is absent in the hole- 
doped side {n < 1). This phase diagram is very similar to 
the one of the three-orbital Hubbard model and also to 
the one of a model for high-Tc cuprates [11] , suggesting 
that the mechanisms of SC are the same origin for both 
systems. It should be noted that the electron-hole asym- 
metry found in a model for high-Tc cuprates, where the 
hole doping favors the SC more than the electron doping, 
is opposite to the results obtained here (see Fig.|4]). This 
difference is simply because of the sign difference of t'/t 
between the effective single-orbital model {t'/t > 0) and 
a model for high-Tc cuprates (t'/t < 0) (STj. 

Although they are qualitatively similar, the phase dia- 
grams for the three-orbital and the effective single-orbital 
Hubbard models are quantitatively different. For in- 
stance, the region of the SC in the single-orbital model 
seems larger than that in the three-orbital model (Fig. 0] 
should be compared with Fig. [2] at U/ti = 8). This in- 
dicates that a multi-orbital effect is to dcstruct the SC. 
The reason can be attributed to the reduction of the ef- 
fective electron correlations due to large orbital fluctua- 
tions, which certainly decrease the probability of facing 
double occupancy at the same orbital. As mentioned 
above, the reduction of the effective electron correlations 
destabilizes the SC and its region becomes smaller. 

In summary, we have studied the ground state phase 
diagram of the three-orbital Hubbard model for Sr2lr04. 
We have found the unconventional SC induced by car- 
rier doping in the Jefi = 1/2 Mott insulator. This SC 
is characterized by the dx2-y2-wavc "pseudospin singlet" 
formed by the Jeff = 1/2 Kramers doublet. We have 
shown that the SC is induced only by electron doping, 
but not by hole doping, for the case of carrier doped 
Sr2lr04. By studying the effective single-orbital Hub- 
bard model constructed from the "Jeff = 1/2 band", we 
have found the similar phase diagram to the one of a 
model for high-Te cuprates, suggesting the same mech- 
anism of the SC in both systems. Finally, it should be 
noted that SC has not been observed yet experimentally 
in layered perovskite Ir oxides. We hope that our study 
will stimulate further experimental as well as theoretical 
studies in this direction. 
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